TITLE OF THE INVENTION 

METHOD OF AND DEVICE FOR SIMULATION 

BACKGROUND OF THE INVENTION 
5 Field of the Invention 

The present invention relates to a method of and device for simulation which 
generates models that reflect statistical variations (the amoimt of change) in electrical 
characteristics of a device for determination of margins in circuit design. 

1 0 Description of the Background Art 

In recent years, the size reduction of devices has been rapidly promoted, and the 
degree of integration of semiconductor integrated circuits has accordingly been 
dramatically on the increase. At the present time, simulation devices such as SPICE 
(Simulation Program with Integrated Circuit Emphasis) are used to assist in designing 

1 5 semiconductor integrated circuits having such an increased degree of integration. 

Conventionally, there are variations in statistics on values (referred to 
hereinafter as "model parameters") regarding the shape of devices in manufactured 
semiconductor integrated circuits and parameters (referred to hereinafter as "process 
parameters"; the model parameters and the process parameters are together referred to as 

20 "device parameters") regarding manufacturing conditions during manufacturing steps 
such as the concentration of impurity for implantation into semiconductor substrates. 
Such variations in model parameters and process parameters give rise to variations in 
electrical characteristics of the devices during operation. With MISFETs (Metal 
Insulator Semiconductor Field Effect Transistors) as an example, variations in model 

25 parameters such as a gate length and a gate insulation film thickness give rise to 



variations in threshold voltage value and in saturation current value in constant current 
regions. 

As the device size reduction is promoted, the variations in model parameters 
and process parameters are not necessarily scaled down (or reduced) in a similar maimer. 

5 Thus, the influence of the variations in these parameters upon the variations in electrical 
characteristics is becoming non-negligible. 

In circuit design with suitable margins, it is therefore important to cause a 
simulation device to reflect the variations in model parameters and process parameters 
and the influence thereof upon the variations in electrical characteristics of the devices. 

10 Although various methods of such reflection have been attempted, attention will now be 
directed to a method which uses a model generally known as a "comer model." 

The comer model refers to a model wherein some of the varying electrical 
characteristics of a device which show a lower j&equency of occurrence than statistical 
typical values (average values, median values and the like) are expressed as comers 

1 5 which define the limits of the variations. At the comers, to what extent the variations in 
model parameters and process parameters are tolerated is calculated. This comer model 
will be described, taking a CMOS (Complementary Metal Oxide Semiconductor) inverter 
as an example. 

Examples of the model parameters used herein include a channel length 
20 difference AL between the geometric gate length and the effective channel length of 
NMOS and PMOS transistors constituting the CMOS inverter, a chaimel width difference 
AW between the geometric gate width and the effective channel width of the transistors, 
and a gate insulation film thickness Tox. The geometric gate length and the geometric 
gate width refer to the gate length and the gate width W^, respectively, of a gate 
25 electrode 5 as viewed in outside configuration in a MISFET structure shown in cross- 
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section of Fig. 7 and in top plan of Fig. 8. The effective channel length and the effective 
channel width refer to the effective length and w^idth, respectively, of a channel formed in 
a channel region 1 1 between source/drain structures 2 and 3. 

An example of the process parameters used herein includes a threshold voltage 
5 VthO at the time that a body voltage is 0 V, the threshold voltage VthO being determined 
by an impurity concentration and the like. Examples of the electrical characteristics of a 
device used herein include a saturation current value Idsat in a constant current region and 
a threshold voltage value Vth. 

It is assumed that the CMOS inverter has the highest response speed in a 

1 0 condition A and has the lowest response speed in a condition B. Considering the NMOS 
and PMOS transistors as a pair, there are four possible conditions: a condition AA 
wherein both the NMOS and PMOS transistors are in the condition A; a condition BB 
wherein both the NMOS and PMOS transistors are in the condition B; a condition AB 
wherein the NMOS and PMOS transistors are in the conditions A and B, respectively; and 

15 a condition BA wherein the NMOS and PMOS transistors are in the conditions B and A, 
respectively. 

Fig. 9 is a graph showing an NMOS saturation current value Idsatn on the 
horizontal axis versus a PMOS saturation current value Idsatp on the vertical axis. In 
the graph of Fig. 9, the point of intersection of the line indicating a typical value Itypn of 

20 the NMOS saturation current value and the line indicating a typical value Itypp of the 
PMOS saturation current value is defined as a point POa. Then, variations in saturation 
current value in the constant current region can be expressed, for example, as a 
quadrilateral SQa surrounding the point POa and having four points PI a to P4a at the 
comers. The area of the quadrilateral SQa and the position thereof relative to the point 

25 POa which fall within a desired range provide semiconductor integrated circuits having 



desired electrical characteristics, to increase the yield of the semiconductor integrated 
circuits. 

Of the four points at the comers of the quadrilateral SQa, the point Pla has an 
NMOS saturation current value greater than the typical value Itypn of the NMOS 

5 satxiration current by (SEdsatna , and a PMOS saturation current value greater than the 
typical value Itypp of the PMOS saturation current by (5Idsatpa . The point P3a has an 
NMOS saturation current value less than the typical value Itypn by (Sldsatnb, and a 
PMOS saturation current value less than the typical value Itypp by Sdsatpb. The 
greater the saturation current value is, the higher the response speed of the NMOS and 

10 PMOS transistors is. Therefore, the point Pla at which both of the saturation current 
values Idsatn and Idsatp are maximized corresponds to the condition AA, and the point 
P3a at which both are minimized corresponds to the condition BB. 

In the light of the steps of forming the CMOS inverter, it is difficult to imagine 
the occurrence of a combination of an NMOS saturation current value greater than the 

1 5 typical value Itypn by tJldsatna and a PMOS saturation current value less than the typical 
value Itypp by fSIdsatpb. Therefore, the point P4a having a smaller difference in 
characteristic value between the NMOS and PMOS transistors than a point P5a indicating 
the above-mentioned combination corresponds to the condition AB. Similarly, the point 
P2a having a smaller difference in characteristic value between the NMOS and PMOS 

20 transistors than a point P6a corresponds to the condition BA, 

The determination of a set of variations from the device parameter values of a 
device corresponding to the point POa (i.e., a set of differences 5(AL) , S{AW) , ^^ox 
and <5VthO between the device parameters values at the point POa and the device 
parameter values at each of the comer points; which are referred to hereinafter as a device 

25 parameter set) at each of the comer points provides the range of the variations in the 
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device parameters corresponding to the quadrilateral SQa indicating the entire area of the 
variations in the saturation current values Idsatn and Idsatp. 

Fig. 10 is a graph showing the absolute value Vthn of an NMOS threshold 
voltage on the horizontal axis versus the absolute value Vthp of a PMOS threshold 

5 voltage on the vertical axis. In the graph of Fig. 10, the point of intersection of the line 
indicating a typical value Vtypn of the absolute value of the NMOS threshold voltage and 
the line indicating a tj^ical value Vtypp of the absolute value of the PMOS threshold 
voltage is defined as a point POb. Then, variations in threshold voltage can be expressed 
as a quadrilateral SQb surrounding the point POb and having four points Plb to P4b at the 

10 comers, in a manner similar to the quadrilateral SQa. 

Of the four points at the comers of the quadrilateral SQb, the point Plb has an 
NMOS threshold voltage value greater than the typical value Vtypn of the absolute value 
of the NMOS threshold voltage by c5Vthna , and a PMOS threshold voltage value greater 
than the typical value Vtypp of the absolute value of the PMOS threshold voltage by 

15 <5Vthpa. The point P3b has an NMOS threshold voltage value less than the typical 
value Vtypn by <5Vthnb , and a PMOS threshold voltage value less than the typical value 
Vtypp by SVthph. 

The greater the absolute value of the threshold voltage is, the lower the 
response speed of the NMOS and PMOS transistors is. Contrary to the points of the 
20 quadrilateral SQa, the point Plb at which both of the absolute values Vthn and Vthp of 
the threshold voltages are maximized corresponds to the condition BB, and the point P3b 
at which both are minimized corresponds to the condition AA. The point P2b 
corresponds to the condition AB, and the point P4b corresponds to the condition B A. 

The determination of a device parameter set at each of the comer points 
25 provides the range of variations in the device parameters corresponding to the 
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quadrilateral SQb indicating the entire area of the variations in the absolute values Vthn 

and Vthp of the threshold voltages. 

Conventionally, a device parameter set at each comer point has been 

determined by a method as shown in Fig. 1 1 . First, initial values of a device parameter 
5 set at each comer point varying depending on variations are provided to a simulator as in 

step SI 02. The initial values used herein may include recorded values of the variations 

specified based on past data in a manufacturing line, values of measurement obtained by 

new sampling, and the like. 

Then, as in Step 8103, desired electrical characteristic values (e.g., Idsatn + 
10 ^adsatna and Vthn + ^thna in the above instances) tolerable for each comer and 

information (a gate length and a gate width; and a body potential and a gate bias value) 

about a device of interest are provided to the simulator. 

Next, using the provided initial values of the device parameter set, the simulator 

performs a circuit simulation based on the information about the device of interest as in 
15 Step SI 04 to calculate electrical characteristic values corresponding to the device 

parameter set. 

Subsequently, as in Step SI 05, whether or not the electrical characteristic 
values satisfy the desired electrical characteristic values provided in Step SI 03 is verified. 
If they do, the device parameter set is adopted as a device parameter set at that comer as 
20 in Step SI 07. 

If they do not, the values of sensitive device parameters are re-adjusted in 
accordance with the device size as in Step SI 06, and the circuit simulation is performed 
again. The operations of simulation, verification and device parameter set re-adjustment 
are repeated until the calculated electrical characteristic values satisfy the tolerable 
25 electrical characteristics provided in Step S 1 03 . 



The expression "sensitive device parameters in accordance with the device size" 
refers to device parameters listed in Table 1, for example, for a MISFET. In Table 1, 
whether the transistor gate size is long or short is determined depending on whether or not 
an electrical characteristic is sensitive to a slight variation in channel length, and whether 
5 the transistor gate size is wide or narrow is determined depending on whether or not the 
electrical characteristic is sensitive to a slight variation in channel width. 



Table 1 



transistor size 


parameter 


LONG/WIDE 


Tox, VthO 


SHORTAVIDE 


AL 


SHORT/NARROW 


AW 



10 The determination of a device parameter set at the comer point P2a in the 

instance of Fig. 9 will be taken as an example for description. First, initial values of the 
device parameter set are provided to a simulation device (Step SI 02). A desired 
electrical characteristic value tolerable for the point F2a is also provided to the simulation 
device (Step 8103), 

15 Then, a circuit sunulation is performed (Step SI 04). It is assiimed that the 

electrical characteristic at a point P7a in Fig. 9 is determined. Subsequently, a 
judgement is made as to whether or not the electrical characteristic value at the point P7a 
falls within a range regarded as approximately equal to the electrical characteristic value 
at the point P2a (Step SI 05). If judged as being within the range, the electrical 

20 characteristic value at the point P7a is adopted as the electrical characteristic value at the 
point P2a (Step SI 07). 
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If the electrical characteristic value at the point P7a is not judged as being 
within the range, the re-adjustment of the sensitive parameter values and the circuit 
simulation are repeated until the electrical characteristic value falls within the range 
regarded as approximately equal to the electrical characteristic value at the point P2a 
5 (StepS106). 

A similar process is performed to sequentially determine device parameter sets 
at the remaining comer points Pla, P3a and P4a. 

Unfortunately, the above-mentioned simulation method, in which the circuit 
simulation calculation, verification and device parameter set re-adjustment are repeated 
10 until the desired electrical characteristic values are reached, requires much time and 
handling to determine the values of the device parameter set at each comer. 

Moreover, which device parameter value is to be adjusted in the device 
parameter set and to what extent it is to be adjusted are not clear. Therefore, the above- 
mentioned simulation method is required to gradually re-adjust the device parameter set 
15 by trial and error, presenting difficulties in determining all of the adjusted values of the 
device parameters at once (or uniquely). 



SUMMARY OF THE INVENTION 

A first aspect of the present invention is intended for a method of simulation, 

20 wherein variations in an electrical characteristic of a device constituting a semiconductor 
integrated circuit are represented in the form of a comer model including at least one 
comer defming a limit of the variations. According to the present invention, the method 
comprises the steps of: (a) preparing a predetermined value tolerable for the variations in 
the electrical characteristic at the at least one comer; (b) performing a circuit simulation 

25 to determine a device parameter sensitivity which is the derivative of the electrical 
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characteristic with respect to a device parameter indicative of information about the 
device; and (c) applying the device parameter sensitivity and the predetermined value of 
the electrical characteristic to the normal equation of the least squares method to 
determine variations in the device parameter at the at least one comer. 
5 Preferably, according to a second aspect of the present invention, the method of 

the first aspect further comprises the step of (d) calculating the variations in the electrical 
characteristic at the at least one comer, based on the multiplication of the device 
parameter sensitivity provided in the step (b) and the variations in the device parameter at 
the at least one comer provided in the step (c). 

10 Preferably, according to a third aspect of the present invention, in the method of 

the second aspect, a comparison is made between the variations in the electrical 
characteristic at the at least one comer calculated in the step (d) and the predetermined 
value prepared in the step (a), and, if an error of the electrical characteristic is greater than 
another predetermined value, the steps (b) through (d) are executed again. 

1 5 Preferably, according to a fourth aspect of the present invention, in the method 

of the second aspect, a comparison is made between the variations in the electrical 
characteristic at the at least one comer calculated in the step (d) and the predetermined 
value prepared in the step (a), and, if an error of the electrical characteristic is greater than 
another predetermined value, a new device parameter is introduced to execute the steps 

20 (b) through (d) using the new device parameter and the device parameter in combination. 
Preferably, according to a fifth aspect of the present invention, in the method of 
the fnst aspect, the device parameter includes a plurality of device parameters; and the 
step (c) is not executed upon at least one of the device parameters, but is executed upon 
only the remainder of the device parameters. 

25 Preferably, according to a sixth aspect of the present invention, in the method of 
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the first aspect, the variations in the device parameter are determined using the weighted 
least squares method in the step (c). 

A seventh aspect of the present invention is intended for a device for simulation, 
the device using a method of simulation as recited in any one of the first to sixth aspects 
5 to represent the variations hi the electrical characteristic of the device in the form of the 
comer model. According to the present invention, the device comprises: data input 
means for inputting the predetermined value tolerable for the variations in the electrical 
characteristic; data output means; a simulator for simulating the amount of change in the 
electrical characteristic as would occur when the device parameter is changed, to 

1 0 determine the device parameter sensitivity; and data processing means for applying the 
device parameter sensitivity determined by the simulator and the predetermined value of 
the electrical characteristic inputted to the data input means to the normal equation of the 
least squares method, to determine variations in the device parameter at the at least one 
comer, thereby outputting the variations in the device parameter at the at least one comer 

15 to the data output means. 

In accordance with the first aspect of the present invention, the variations in the 
device parameter at the at least one comer are determined by applying the device 
parameter sensitivity and the predetermined value of the electrical characteristic to the 
normal equation of the least squares method. Therefore, the method can determine the 

20 variations in the device parameter tolerable for the at least one comer without the need to 
repeat the circuit simulation and also can uniquely determine the variations in the device 
parameter. 

In accordance with the second aspect of the present invention, the variations in 
the electrical characteristic at the at least one comer are calculated based on the 
25 multiplication of the device parameter sensitivity and the variations in the device 
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parameter at the at least one comer. Therefore, the method can verify whether or not the 
value of the variations in the electrical characteristic calculated in the step (d) is regarded 
as approximately equal to the predetermined value prepared in the step (a). 

In accordance with the third aspect of the present invention, the comparison is 

5 made between the variations in the electrical characteristic at the at least one comer 
calculated in the step (d) and the predetermined value prepared in the step (a), and, if the 
error of the electrical characteristic is greater than another predetermined value, the steps 
(b) through (d) are executed again. Therefore, the method might determine the device 
parameter and the electrical characteristic at the at least one comer with high accuracy. 

10 In accordance with the fourth aspect of the present invention, the new device 

parameter is introduced, and the steps (b) through (d) are executed using the new device 
parameter and the device parameter in combination. This further increases the accuracy 
of the electrical characteristic provided in the step (d). 

In accordance with the fifth aspect of the present invention, the step (c) is not 

15 executed upon at least one of the device parameters but is executed upon only the 
remainder of the device parameters. Therefore, the method can determine the variations 
in the device parameter tolerable for the at least one comer while reducing the amount of 
calculation. 

In accordance with the sixth aspect of the present invention, the weighted least 
20 squares method is used to determine the variations in the device parameter in the step (c). 
Therefore, the method can determine the variations in the device parameter more 
correctly while emphasizing the error of the electrical characteristic to which particular 
importance is desired to be attached. 

In accordance with the seventh aspect of the present invention, there is provided 
25 a device for simulation which implements the method of simulation as set forth in any one 



12 

of the first to sixth aspects of the present invention. 

It is therefore an object of the present invention to provide a method of and 
device for simulation which can determine the values of a device parameter set at each 
comer without repeating a circuit simulation and which can xmiquely determine the values 
5 of the device parameter set. 

These and other objects, features, aspects and advantages of the present 
invention will become more apparent from the following detailed description of the 
present invention when taken in conjunction with the accompanying drawings. 

1 0 BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a flowchart showing a method of simulation according to a first 
preferred embodiment of the present invention; 

Fig. 2 is a flowchart showing the method of simulation according to a second 
preferred embodiment of the present invention; 
15 Fig. 3 is a flowchart showing the method of simulation according to a third 

preferred embodiment of the present invention; 

Fig. 4 shows a simulation device according to a fourth preferred embodiment of 
the present invention; 

Figs. 5 and 6 are flowcharts showing the method of simulation according to the 
20 fourth preferred embodiment; 

Fig. 7 is a cross-sectional view of a MISFET structure; 
Fig. 8 is a top plan view of the MISFET structure; 
Figs. 9 and 10 illustrate comer models; and 

Fig, 1 1 is a flowchart showing a background art method of simulation. 

25 
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DESCRIPTION OF THE PREFERRED EMBODIMENTS 
<First Preferred Embodiment> 

A first preferred embodiment of the present invention is intended for providing 
a method of simulation w^hich can determine the values of a device parameter set at each 
5 comer without repeating a circuit simulation and which can uniquely determine the values 
of the device parameter set. To provide the method, the present invention uses the linear 
least squares method. 

The first preferred embodiment will now be described based on the comer 
models for a CMOS inverter shown in Figs. 9 and 10 as an example, as described in the 

10 background art. It is assumed that m CMOS inverters are present in the first preferred 
embodiment. As in the background art, examples of the model parameters used herein 
include a channel length difference AL between the geometric gate length and the 
effective channel length of NMOS and PMOS transistors constituting the CMOS inverters, 
a channel width difference AW between the geometric gate width and the effective 

15 channel width of the transistors, and a gate insulation film thickness Tox. An example 
of the process parameters used herein includes a threshold voltage VthO at the time that a 
body voltage is 0 V. The threshold voltage VthO at the time the body voltage is 0 V may 
be, for example, that of a transistor having a gate size which is "LONG/WIDE." 
Examples of the electrical characteristics of the devices used herein include a saturation 

20 current value Idsat in a constant current region and a threshold voltage value Vth. 

In regions (the region indicated by the quadrilateral SQa of Fig. 9 and the 
region indicated by the quadrilateral SQb of Fig. 10) wherein comer models are to be 
produced, variations in electrical characteristics may be regarded as linear with variations 
in device parameter set. The variations in electrical characteristics near the points POa 

25 and POb and the corresponding variations in device parameter set are slight enough to be 
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regarded as changing linearly. 

The relationships between the variations in electrical characteristics and the 
variations in device parameter set for the NMOS and PMOS transistors are expressed by 
Equations (1) through (4) in which data about the electrical characteristics to be observed 
by simulation is on the left-hand side and terms of device parameters responsible for 
changes in the electrical characteristics are on the right-hand side. 



9Idsatn ^ SIdsatn 

Itypn + adsatn = Itypn + ALn) — — — + S ( AWn) 

^ SIdsatn ^ SIdsatn ^ 

+ (JToxn + (5VthOn— — — + fin 

aToxn aVthOn 



(1) 



aVthn ,5Vthn 
Vtypn + ^thn = Vtypn + J(ALn) + ^(AWn) 



aALn SAWn 



5Vthn „ , , ^ 9Vthn „^ 

+ SToxn + <5VthOn -— — + f 2n 

SToxn SVthOn 



^ _ ^ , aidsatp A N SIdsatp 
Itypp + 5[dsatp = Itypp + ^(ALp) + J(AWp) 

^ SIdsatp ^ SIdsatp , 

+ t^oxp + tWthOp — — r + f 

^ aroxp avthop 



... (3) 



9VthD 9Vthp 
Vtypp + ^thp - Vtypp + ^( ALp) — — ^ + S ( AWp) — 

9ALp 5AWp 

^ aioxp avthOp 



where Itypn is a typical value of an NMOS saturation current value at the point POa, Itypp 
is a typical value of a PMOS saturation current value at the point POa, Vtypn is a typical 
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value of an NMOS threshold voltage at the point POb, and Vtypp is a typical value of a 
PMOS threshold voltage at the point POb. 

In Equation (1), «5[dsatn , like ^dsatna and <5[dsatnb shown in Fig. 9, is a 
variation from the typical value Itypn; <5(ALn) is a variation from a chaimel length 
5 difference ALn in the device parameter set corresponding to the NMOS electrical 
characteristics at the point POa; and ^(AWn) , ^oxa and JVthOn are variations 
defined in a like manner. Also, Sldsatn/aALn , Sldsatn/SAWn, aidsatn/SAToxn 
and aidsatn/aAVthOn which are the derivatives of the electrical characteristics with 
respect to the device parameters are device parameter sensitivities of the NMOS 
1 0 saturation current value to the device parameters, respectively. 

Likewise, the terms in Equations (2) through (4) denote variations in the 
saturation current value and the threshold voltage from their typical values for the NMOS 
and PMOS transistors, variations in device parameters corresponding to the typical values, 
and the device parameter sensitivities of the saturation current value and the threshold 
15 voltage to the device parameters. The terms fin, f2n, flp, f2p on the right-hand side of 
Equations (1) through (4) are terms other than those regarding the device parameters and 
denoting errors such as calculation errors during simulation. 

A vector x, a vector 0 , a matrix A and a vector f are defined as 
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x = 



(Sdsatnl ^ 

Sdsatnm 
<5Vthnl 

<5Vtlmm 
Sdsatpl 

<5Idsatpm 
(5Vthpl 

^thpm 



... (5) 



e = 



f SiALn) ^ 
<5(AWn) 
JToxn 
^thOn 
5(ALp) 
<5(AWp) 
JToxp 
(5VthOp 



... (6) 



eidsatnl ai dsatnl aidsatnl gldsatnl 

aVthOn 

aVthnl 



0 



0 



0 



0 



5ALn 9AWn . eioxn 
avthnl aVthnl ' SVthnl 



5ALn 



0 



aAWn 
0 

0 



SToxn 
0 

0 



aVthOn 
0 

0 



0 0 0 0 

aidsatpl aidsatpl aidsatpl aidsatpl 

aALp aAWp aioxp avthOp 

avthpi avthpi " avthpi avthpi 
aALp aAWp aioxp avthOp 

... (7) 
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r f Inl ^ 

flrnn 
f2nl 



f = 



f2nm 
flpl 



flpm 
f2pl 



(8) 



Then, Equations (1) through (4) axe expressed together as 

x = A<9 + f -(9) 

wherein Itypn, Vtypn, Itypp and Vtypp which are present on both sides are eliminated. 

Because of the presence of m CMOS inverters, m equations are formulated for each of 

Equations (1) through (4). This is indicated by the adscripts 1 to m added to Idsatn, 

Vthn, fin, f2n and the like in the vector x, the matrix A and the vector f. 

According to the linear least squares method, an observed value error is 

minimized when the square of the norm of the vector f of an error has the least value. 

Therefore, it is required to determine the inflection point of the value of the square of the 

norm of 

f =^(yspeci-ysimi)^ 

i=l 
N 

= ^ (^speci - (^'simi)^ 

i=l 

='(x-A^)(x-A^) 

... (10) 
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where yspeci in the first expression on the right-hand side is an electrical characteristic 
value at a comer point tolerable for variations, ysimi is an electrical characteristic value at 
a comer point observed by simulation, and i = 1 to N (N = 4 herein). For example, 
yspeci is obtained by subtracting fin and Itypn from the right-hand side of Equation (1); 
5 yspec2 is obtained by subtracting fZn and Vtypn from the right-hand side of Equation (2); 
yspec3 is obtained by subtracting flp and Itypp from the right-hand side of Equation (3); 
yspec4 is obtained by subtracting f2p and Vtypp from the right-hand side of Equation (4); 
ysiml is obtained by subtracting Itypn from the left-hand side of Equation (1); ysim2 is 
obtained by subtracting Vtypn from the left-hand side of Equation (2); ysimB is obtained 
10 by subtracting Itypp from the left-hand side of Equation (3); and ysim4 is obtained by 
subtracting Vtypp from the left-hand side of Equation (4). 

In Equation (10), ^speci in the second expression on the right-hand side 
corresponds to components of the product of the matrix A and the vector d in the first 
term on the right-hand side of Equation (9); <%^simi corresponds to components of the 
1 5 vector x on the left-hand side of Equation (9); and ' (x - A^) in the third expression on 
the right-hand side is the transpose of the vector (x - A^) . 

The inflection point of the value of the square of the norm of Equation (10) is 
determined by solving 



20 The partial differentiation with respect to the vector d in Equation (11) means the 
differentiation with respect to each of the components of the vector 0 . Since the vector 
0 used herein has eight components, Equation (11) means substantially eight 
simultaneous equations. 

Using the third expression on the right-hand side of Equation (10), Equation 
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(1 1) is solved, with the matrix representation maintained. 

^AA^=*Ax ".(12) 

Equation (12) is referred to as the normal equation of the least squares method. In 

Equation (12), *A is the transpose of the matrix A. 

5 If *AA on the left-hand side of Equation (12) is regular, the least squares 

estimate is given by 

^ = (^Aa)'' ^Ax ...(13) 

where the vector ^e has components which are the estimates of the respective 
components of the vector 0 . 
10 Thus, determining the device parameter sensitivities which are the components 

of the matrix A and calculating the vector at each comer point using Equation (13) 
provide a device parameter set corresponding to the electrical characteristics required for 
each comer point. 

The device parameter sensitivities are determined by performing a circuit 
15 simulation. More specifically, the device parameter sensitivities are determined by 
observing the amount of slight change in electrical characteristics from the electrical 
characteristics at the points POa and POb which are provided when the components of the 
device parameter set are slightly changed one by one. 

The above-mentioned procedure will be described with reference to Fig. 1. 
20 Fig. 1 is a flowchart showing the method of simulation according to the first preferred 
embodiment of the present invention. 

First, values (Itypn + ^dsatna , Itypn - <5[dsatnb and the like) tolerable for 
variations in electrical characteristics at each comer and information (the typical values of 
the model and process parameters corresponding to the points POa and POb) about a 
25 device of interest are prepared (Step S02). Then, a circuit simulation is performed to 
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determine the device parameter sensitivities of the electrical characteristics to the device 

parameters (Step S03). 

Equation (13) derived from the normal equation of the least squares method is 

used to calculate variations iS(AL) , ^(AW) , ^ox , SVihO ) in device parameter set at 
5 each comer from the device parameter sensitivities provided in Step SOS and the electrical 

characteristic values provided in Step S02 (Step S04). 

The simulation method according to the first preferred embodiment determines 

the variations in device parameter set at each comer by applying the device parameter 

sensitivities and the electrical characteristic values to the normal equation of the least 
10 squares method, to require the circuit simulation to be performed only when determining 

the device parameter sensitivities. This eliminates the need to repeat the circuit 

simulation to determine the variations in the device parameters tolerable for the comers. 

Additionally, Equation (11) which means a system of simultaneous equations the number 

of which equals the number of device parameters is uniquely solved based on a 
15 fundamental theorem of algebra. Therefore, the simulation method of the first preferred 

embodiment can uniquely determine the variations in the device parameters. 

<Second Preferred Embodiment> 

A method of simulation according to a second preferred embodiment of the 
20 present invention is a modification of the method of simulation according to the first 
preferred embodiment. In the method according to the second preferred embodiment, 
the values of variations in device parameter set at each comer which are determined in the 
first preferred embodiment are used to actually calculate the electrical characteristics at 
each comer and to verify whether or not the calculated electrical characteristics are 
25 regarded as approximately equal to the established values tolerable for electrical 
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characteristic variations. 

Fig. 2 is a flowchart showing the method of simulation according to the second 
preferred embodiment. First, as in the first preferred embodiment, the values tolerable 
for variations in electrical characteristics at each comer and information about a device of 
5 interest are prepared (Step S12). 

Then, a circuit simulation is performed to determine the device parameter 
sensitivities of the electrical characteristics to the device parameters (Step SI 3). 
Equation (13) derived from the normal equation of the least squares method is used to 
calculate variations in device parameter set at each comer from the device parameter 
10 sensitivities provided in Step S 13 and the electrical characteristic values provided in Step 
S12(StepS14). 

Thereafter, the electrical characteristic values at each comer are calculated 

based on the variations in device parameter set provided in Step S14 (Step SI 5). 

Equation (9) from which the vector f is eliminated should be used for the calculation of 
15 the electrical characteristic values. In other words, the variations in electrical 

characteristics should be calculated based on the multiplication of the device parameter 

sensitivities and the device parameter set. 

Then, whether or not the electrical characteristic values provided in Step SI 5 

satisfy the established values tolerable for electrical characteristic variations is verified 
20 (Step S 16). This verification may be performed, for example, using the value of an error 

ERR determined by 

2 _ (ysimin-yspecin)- +(ysimip-yspecip)^ ^^^^ 
yspecin^ +yspecip" 
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ERR = .^^CERRff' + ERRss' + ERR^s' + ERR^p' ) ... (15) 

where ERRxx^ is the square error of the electrical characteristics at each comer (XX is a 
subscript indicative of NMOS and PMOS response speeds, respectively; F indicates the 
fastest response speed and S indicates the slowest response speed), and ysimin, ysimip, 
5 yspecin, yspecip are the values of ysimi and yspeci in the first expression on the right- 
hand side of Equation (10) for NMOS and PMOS transistors. 

If the value of the error ERR is greater than a predetermined value, there is a 
possibility that, for example, the calculation of the device parameter sensitivities has an 
error. Then, the flow should be returned to Step S13 wherein the device parameter 
10 sensitivities are calculated again. This might reduce the error of the device parameter 
sensitivities to make the value of the error ERR less than the predetermined value, in 
which case the device parameter set and the electrical characteristics with high accuracy 
are provided. 

The simulation method according to the second preferred embodiment 
15 calculates the variations in the electrical characteristics at each comer based on the 
multipUcation of the device parameter sensitivities and the variations in device parameters 
at each comer, to verify whether or not the calculated electrical characteristics are 
regarded as approximately equal to the established values tolerable for electrical 
characteristic variations. 

20 Additionally, the flow returns to Step S13 to calculate the device parameter 

sensitivities again if the value of the error ERR is greater than the predetermined value. 
This might make the value of the error ERR less than the predetermined value to provide 
the device parameter set and the electrical characteristics with high accuracy. 
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<Third Preferred Embodiment> 

A method of simulation according to a third preferred embodiment of the 
present invention is a modification of the method of simulation according to the second 
preferred embodiment. The third preferred embodiment differs from the second 
preferred embodiment in the step which follows the step of verifying whether or not the 
calculated electrical characteristics at each comer are regarded as approximately equal to 
the established values tolerable for the electrical characteristic variations. If the 
calculated electrical characteristics are not regarded as approximately equal to the 
established values, the dimensions of the matrix A and the vector d in Equation (9) are 
adjusted or the device parameters are weighted. 

Fig. 3 is a flowchart showing the method of simulation according to the third 
preferred embodiment. First, as in the second preferred embodiment, the values 
tolerable for variations in electrical characteristics at each corner and information about a 
device of interest are prepared (Step S22). 

Then, a circuit simulation is performed to determine the device parameter 
sensitivities of the electrical characteristics to the device parameters (Step S23). 
Equation (13) derived from the normal equation of the least squares method is used to 
calculate variations in device parameter set at each comer from the device parameter 
sensitivities provided in Step S23 and the electrical characteristic values provided in Step 
S22 (Step S24). 

Thereafter, the electrical characteristic values at each comer are calculated 
based on the variations in device parameter set provided in Step 824 (Step 825). Then, 
whether or not the electrical characteristic values provided in Step 825 satisfy the 
estabhshed values tolerable for electrical characteristic variations is verified (Step S26), 
This verification may be performed, for example, using the value of the error ERR, as in 
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the second preferred embodiment. 

If the value of the error ERR is greater than a predetermined value, the 
dimensions of the matrix A and the vector d in Equation (9) are adjusted or the device 
parameters are weighted. 

The expression "adjusting the dimensions of the matrix A and the vector d in 
Equation (9)" refers to adding a new device parameter which influences the electrical 
characteristic values to increase the dimension of the matrix A and the vector 6 . In 
this case, the addition of the new device parameter involves the need to newly calculate 
the device parameter sensitivities of the electrical characteristics to the new device 
parameter in Step S23 (which is indicated by the arrow directed from Step S27 to Step 
S23 in Fig. 3). 

That is, the new device parameter is introduced, and the calculation of the 
device parameter sensitivities and the application of the least squares method are carried 
out upon the new device parameter. This increases the accuracy of the electrical 
characteristics provided in Step S25, 

The expression "weighting the device parameters*' refers to modifying Equation 

(13) into 

^ = (^AWa)~' *AWx ...(16) 
where W is a weighting matrix defined by 







0 


0 


0 ^ 




0 


W2 


0 • 


• 0 


w = 


0 


0 


W3 


0 




.0 


0 


0 





The weighting matrix W is a diagonal matrix having diagonal components Wl to Wn to 
be assigned as weights to the components of the vector x and the components in the 
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respective rows of the matrix A. This process corresponds to multiplying a weight by 
^speci and (5ysimi in the second expression on the right-hand side of Equation (10) for 
each electrical characteristic, and also means changing the degree of contribution of each 
electrical characteristic to the norm of the vector f. Such a technique is known as a 
5 weighted least squares method which can determine the device parameter set more 
correctly while emphasizing the error of the component of an electrical characteristic to 
which particular importance is desired to be attached among the components of the vector 

X. 

Although the weighting process is performed in Step 827 after the step of 
10 verification in Step S26 in the third preferred embodiment, the device parameters may be 
weighted in Step S24 when the device parameter set is calculated. Altematively, even if 
the calculated electrical characteristics at each comer are regarded as approximately equal 
to the established values tolerable for the electrical characteristic variations in Step S26, 
the device parameters may be weighted and then Step S24 may be executed (which is 
15 indicated by the broken arrow directed from the arrow with YES in Step S26 to Step 
S27). 

Further, even if the calculated electrical characteristics at each comer are 
regarded as approximately equal to the established values tolerable for the electrical 
characteristic variations in Step S26, an adjustment may be made such that a device 

20 parameter of a low parameter sensitivity, for example, is eliminated from the matrix A 
and the vector 9 in Equation (9) in Step S27 to decrease the dimensions of the matrix A 
and the vector 9 (e.g., the component ^oxn in the vector Q is removed and the 
components in the third column of the matrix A are removed), and only the device 
parameters of a high device parameter sensitivity are used to calculate the device 

25 parameter set again (which is also indicated by the broken arrow directed from the arrow 



with YES in Step S26 to Step S27). In other words, without applying one or some of the 
device parameters to the least squares method. Step S24 is performed again only on the 
remaining device parameters. This reduces the amount of calculation to rapidly 
calculate the device parameter set and the electrical characteristics again in Steps S24 and 

5 S25 while placing stress on the device parameters of a high device parameter sensitivity. 

The simulation method according to the third preferred embodiment adjusts the 
dimensions of the matrix A and the vector 6 in Equation (9) to further increase the 
accuracy of the electrical characteristics. Otherwise, the method can reduce the amount 
of calculation to rapidly calculate the device parameter set in Step S24 while placing 

10 stress upon the device parameters of a high device parameter sensitivity. Additionally, 
when the weighted least squares method is used, the method of the third preferred 
embodiment can determine the device parameter set more correctly while emphasizing 
the error of an electrical characteristic to which particular importance is desired to be 
attached. 

15 

<Fourth Preferred Embodiment> 

A fourth preferred embodiment of the present invention presents a simulation 
device which implements the method of simulation, for example, according to the third 
preferred embodiment. 

20 Fig. 4 shows a construction of the simulation device according to the fourth 

preferred embodiment. The simulation device comprises a data input section 100 for 
inputting numerical data and the like, a data processing section 200 for calculation and 
transfer of various data and the like, a data output section 300 for displaying the results of 
calculation of the device parameter set and the electrical characteristics and other data, a 

25 data storage section 400 for storing various data therein, and a simulator 500 for 
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performing a circuit simulation. 

The data processing section 200 includes, for example, a CPU (Central 
Processing Unit), and comprises an internal storage section 201 such as a cache memory. 
The data storage section 400 is an external storage device, such as a hard disk, which 
5 stores therein a sensitivity file 401 in which the device parameter sensitivities calculated 
in the past are recorded, an extraction result file 402 in which the device parameter sets 
calculated in the past are recorded, an extracted condition file 403 in which extraction 
conditions (MISFET gate bias value, body voltage value and the like) used in calculation 
of the device parameter sets are recorded, and a file 404 in which typical parameter sets 
10 which are sets of device parameters having typical values at the respective points POa, 
POb and the like are recorded. 

Figs. 5 and 6 are flowcharts in greater detail when the flowchart of Fig. 3 is 
applied to the simulation device of the fourth preferred embodiment. 

First, in Step S52, the typical parameter sets, the extraction conditions, 
15 transistor sizes such as the gate length and the gate width, and the values tolerable for the 
variations in electrical characteristics at the comers are inputted through the data input 
section 100. The input data for use in calculation of the device parameter sets and the 
electrical characteristics are recorded in the internal storage section 201 or the data 
storage section 400 by the data processing section 200. 
20 Next, the data processing section 200 examines whether device parameter 

sensitivities have been calculated in the past under the same conditions as the input data 
with reference to the various files 401 to 404 in the data storage section 400 (Step S53). 
If the calculation has been made in the past, the record thereof is loaded from the 
sensitivity file 401 and displayed (Step S54). The data processing section 200 uses the 
25 recorded contents to calculate the device parameter sets and the electrical characteristics. 
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If no device parameter sensitivities have been calculated under the same 
conditions in the past, the data processing section 200 provides the input data to the 
simulator 500. Then, the simulator 500 performs a circuit simulation to calculate the 
device parameter sensitivities of the electrical characteristics (Step S55). The result of 

5 calculation is transferred to the data processing section 200. The data processing section 
200 stores the resuU of calculation in the sensitivity file 401 while displaying the result of 
calculation on the data output section 300. 

Next, the data processing section 200 examines whether or not device 
parameter sets have been calculated in the past under the same conditions as the above- 

10 mentioned input data by using the device parameter sensitivities provided in Step S54 or 
S55 with reference to the various files 401 to 404 in the data storage section 400 (Step 
S56). If the calculation has been made in the past, the record thereof is loaded from the 
extraction result file 402 and displayed (Step S57). The data processing section 200 uses 
the recorded contents to calculate the electrical characteristics. 

15 If no device parameter sets have been calculated under the same conditions in 

the past, the data processing section 200 calculates the device parameter sets using 
Equation (13) (Step S58). The data processing section 200 displays the result of 
calculation on the data output section 300. 

Next, the data processing section 200 applies the device parameter sensitivities 

20 and the device parameter sets to Equation (9) from which the vector f is eliminated, to 
calculate the electrical characteristics at the respective comers (Step S59). Then, the 
data processing section 200 verifies whether or not the result of calculation satisfies the 
values tolerable for the variations in the electrical characteristics inputted in Step S52 
(Step S60). 

25 If the result does not satisfy the tolerable values as a result of the verification, 
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the dimensions of the matrix A and the vector 0 in Equation (9) are adjusted (Step S63) 
or the device parameters are weighted (Step S62). After the adjustment of the 
dimensions in Step S63, the flow retums to Step S53 for calculation of new device 
parameter sensitivities. After the weighting of the device parameters in Step S625 the 
5 flow retums to Step S58 for calculation of new device parameter sets. 

The flow indicated by the broken arrow directed from Step S26 to Step S27 in 
Fig. 3 is not shown in Figs. 5 and 6 for the purposes of avoiding complication. It is, 
however, needless to say that such a flow is provided. 

If the result of calculation satisfies the tolerable values as a result of the 
10 verification, the calculated device parameter sets are stored in the extraction result file 
402 (Step S61). Then, the simulation is completed (Step S64). 

The fourth preferred embodiment has been described using the implementation 
of the simulation method according to the third preferred embodiment as an example. 
To implement the simulation method according to the second preferred embodiment, 
15 Steps S62 and S63 in Fig. 6 are omitted, and the flow indicated by the arrow with NO in 
Step S60 is returned to Step S53, To implement the simulation method according to the 
first preferred embodiment, Steps S59, S60, S62 and S63 in Fig. 6 are omitted, and Step 
S61 is executed after Steps S57 and S58. 

While the invention has been described in detail, the foregoing description is in 
20 all aspects illustrative and not restrictive. It is understood that numerous other 
modifications and variations can be devised without departing firom the scope of the 
invention. 



